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GEOS Satellite Tracking Corrections 
for Refraction in the Ionosphere 

John H. Berbert 
Horace C. Parker 
ABSTRACT 

This document compares the analytic formulations at different elevation angles 
and at a frequency of 2-GHz for the ionospheric refraction corrections used on the 
GEOS satellite tracking data. The formulas and ray-trace results for elevations 
greater than 10°, where most satellite tracking is done, differ in AE, AR, and Ar by 
less than 0,4 miliidegrees (1.4 arc-seconds), 12 meters, and 12 cm/sec, respective- 
ly. In comparison to most operational requirements, this is insignificant. However, 
for the GEOS Observation Systems Intercomparison Investigation, these differences 
are equivalent in size to observed differences in system biases for some of the best 
electronic geodetic tracking systems and are probably contributing to the observed 
biases. 

For all the GEOS ionospheric refraction correction formulas examined in this 
report, the range rate refraction correction for AR is the time derivative of the 
formula for AR, as might be expected. It is much easier to derive AR from AR in 
this manner , rather than independently from basic principles as was usually done in 
the earlier documents. 

The ray-trace results and most of the more detailed analytic correction formulas 
show that the ionospheric refraction correction for range rate on an overhead pass is 
a maximum for elevation angles between 15° and 30° and falls off rapidly for both 
higher and lower elevation angles, contrary to the effect of the troposphere and to 
some reports in the literature. 
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measured elevation 
elevation corrected for refraction 
measured range 
range corrected for refraction 
measured range rate 
range rate corrected for refraction 
rate of change of elevation (refer to Appendix B) 
radio frequency 
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height of the satellite = *^x 10 meters 
height of maximum electron density 

ionospheric scale height for the Chapman model ionosphere 
ionospheric scale height for the DC program 
electron density in electrons/m 
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total electron content in a vertical column of area 1 m , N T 
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electron density in electrons/m at h 
refractivity at height h. 

refractivity at height h^, _(N. m < 0 for the ionosphere) 

refractive index 

Napierian base - 2.71828183 

N. He 1 
im 

elevation correction = E - E 

c 

range correction = R - R 

• c • 

range rate correction = R - R 

c 

earth radius = 6,378,166 meters 

ionospheric scale height factor for GPRO arising from parabolic model 
ionosphere. Nominally GPRO uses H =50 km 

thickness parameter for F layer biparabolic profile in Bent model 
ionosphere 

k = 1/H, = decay constant for topside exponential profile in Bent model 
ionosphere 
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R, 


Rh " h m - y^ in Bent model ionosphere 



GEOS SATELLITE TRACKING CORRECTIONS 
FOR REFRACTION IN THE IONOSPHERE 


INTRODUCTION 

The refraction correction formulas for correcting the elevation (E), range 
(R), and range rate (R) tracking data from the GEOS geodetic tracking systems are not 
standardized. Different refraction correction formulas are applied to the same data by 
different users. The different formulas arise from the use of different mathematical 
models of the atmosphere or from different approximations to the same model. 

This document presents some of the different sets of ionospheric re- 
fraction correction formulas now being used on GEOS data and compares the magni- 
tude, at a frequency of 2 GHz, of the corrections to E, R, andR for different ele- 
vation angles. The 2-GHz frequency is close to the frequencies used by the GRARR, 
the Apollo USB, and various other S-band radar tracking systems. This is a continu- 
ation of a study of the comparison of atmospheric refraction correction formulas for 
GEOS data, in support of the GEOS Observation Systems Intercomparison Investigation. 
The tropospheric refraction correction comparisons are given in reference 1, Not 
all of the tracking systems represented in reference 1 are included here because the 
ionospheric refraction correction is not always applied as an analytic correction 
formula. In particular, the Army Secor ranging system and the Navy TRANET Dop- 
pler system use multiple frequency refraction correction techniques. For C-band 
radar tracking data, at a frequency of 5.7 GHz, the ionospheric correction is very 

small and is therefore often neglected or else lumped together with tropospheric pa- 

o • 

rameters. For example, for C-band at 10 , AR = 0.99 meters, AR = 0.13 cm/sec, 
and AE = 0.076 millidegrees (0.27 seconds of arc). 

The sets of refraction correction equations for AE, AR,and AR being 
used by the Goddard Space Flight Center (GSFC) on GEOS data are compared with each 
other in this document. The formulas for AE and AR are also compared with the re- 
sults of ray-tracings through a modified Chapman ionosphere using the Eastern Test 

Range (ETR) REEK program. In addition, the AR refraction correction formulas are 

d (Ar) 

compared with values of — obtained by differentiating the polynomial fitted to the 
curve of the REEK Ar vs time values. Results of ray-tracing through the Bent model 
ionosphere (reference 23) are plotted for comparison. 

In Appendix A, the ionospheric refraction correction equations are first 
given as they appear in the referenced source documents and then are converted to a 
common notation and functional form to simplify the comparison of these equations. 
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CORRECTION FORMULAS FOR AE, AR, AND AR 

Table 1 lists the GSFC GEOS-2 tracking data processing programs and in- 
dicates which programs apply AE, AR, and AR ionospheric refraction corrections. 

Table 1 

GSFC GEOS-2 Tracking Data Processing Programs 


Processing Programs 

AE 

AR 

AR 

Differential Correction, DC, used with GSFC operational 
Orbit Determination Programs. (reference 2) 

Yes 

Yes 

Yes 

Dr. J. J. Freeman Formulations (reference 5) 

Yes 

Yes 

Yes 

GEOS Preprocessor Program, GPRO, used with the 
Geodetic Data Adjustment Program, GDAP, (ref- 
erence 4) and with the Network Adjustment Programs, 
NAP-1 and NAP-2 (reference 19) 

No 

Yes 

Yes 

NONAME, all versions^ prior to MONSTER (8/20/70) which 
corrects for AR and AR using a table look-up for two 
GRARR VHF stations 

No 

No 

No 

GEOS-2 GRARR Validation Program, GEOVAP (ref- 
erence 26) 

No 

Yes 

Yes 

Definitive Orbit Determination System, DODS (ref- 
erence 24). Same formulation as DC Program, ref- 
erence 2, but implementation differs (reference 25) 

Yes 

Yes 

Yes 

NAP-3 (reference 23) 

j 

Yes 

Yes 

Yes 

Goddard Trajectory Determination System, GTDS, 
Techniques not yet finalized 

Yes 

Yes 

Yes 


The analytic corrections in these programs assume that the proper 
correction to the observations E, R, and R for the net effect of the ionosphere is to 
decrease E by the positive quantity AE, to decrease R by the positive quantity AR, 

• • SO 

and to decrease R by the quantity AR. Both R and AR may be positive or negative, 
but both always have the same sign. Thus, the correction AR has the effect of al- 
ways reducing the magnitude of R. To be compatible with these conventions for the 
analytic corrections, the group option in the REEK ray-trace program was always 
used. A discussion of the meaning of the group option in REEK is given in reference 
28, and an analysis of the ionospheric correction for the GRARR system defining 
when to use group and when to use phase corrections, is given in reference 29> The 
Ae, AR, and AR ionospheric refraction correction formulas used with GEOS-1 and 
GEOS-2 data are listed in Tables 2, 3, and 4, respectively. 
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Table 2 

Elevation Refraction Correction Formulas 



















Table 3 

Range Refraction Correction Formulas 



AR. = -N- He esc E = - M n esc E 






























Table 4 

Range Rate Refraction Correction Formulas 



under development 






















The DC, DODS, FREEMAN, and GEOVAP formulas listed in Table 1 
employ as a factor nominal or zero-order correction functions for Ae , Ar, and Ar, 
which are defined as 

-N. H e 1 -M 

E - E = AE. = — ctn E = -r-2. ctn E (1) 

c 10 h t h t 

R - R = AR. = -N. H e 1 esc E = -M esc E (2) 

c 10 im o 

and the derivative of AR. 

10 

R - R = AR.. = +N. H e 1 E ctn E esc E = M E ctn E esc E (3) 
c io im o 

where M = N. He 1 . Note that M is a multiplicative factor in all these zero-order 
o ini o i C h. 

correction functions. It can be shown that M = N. H e 1 \ n t N. dh, when using 

o im © i 

the Chapman model ionosphere. 

These AE. q and AR^ corrections result from the following assumptions: 

(1) Modified Chapman ionosphere 

(2) Flat layered, horizontally homogeneous ionosphere 

(3) Satellite height, h^, large compared to scale height, H 

(4) Negligible ray path bending 

The GPRO and NAP-3 formulas are different in appearance from the DC, 
DODS, Freeman, and GEOVAP formulas due to the choice of a parabolic ionospheric 
model in GPRO by Willmann (references 4 and 22) and a joined topside exponential, 
and bottomside biparabolic, ionospheric model in NAP-3 by Bent (reference 23) com- 
pared to the Chapman model used by DC, DODS, Freeman, and GEOVAP. The GPRO 
and NAP-3 models enable the correction formulas, which involve integration along the 
ray path, to be integrated in closed form, but at the resulting expense of more com- 
plicated formulas . 
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In order to compare the analytic refraction correction equations, we need 

typical values for N. , h . and H. These are obtained as follows: 
im m 

The refraction index, to first order approximation, is given by (reference 11, eq 2,78) 

80.5 N 40.25 N 

therefore, the refractivity at any height is 


N. = /x. 
i i 


40.25 N 



and the refractivity at the height of maximum electron density is 


40.25 N 


N. = u. - 1 = - ■ 
im r im 


em 


(5) 


which provides one of the required multiplicative factors if N and f are known. The 
scale height in kilometers is determined in this study from 


H = | [30 + 0.2 (h m - 200)] 


(6a) 


which provides the other required multiplicative factor if h^ is known. 

The origin and validity of this type of linear relationship between H and 
h is unknown to us, but the first mention we found of it was by Chamow (1959), ref- 
erence 27, who quotes the scale height as 

H = | j30 +0.2 (h m - 200)] (6b) 

This formula for scale height was also quoted in reference 2 (DC program), 
reference 24 (DODS program), and reference 25 (DODS program). 

A similar linear relationship was proposed by Wright (1960, 1961), ref- 
erences 30 and 31, based on ionograms from 11 stations collected during the period 
January 1959 to February 1960. 

Freeman (1965), reference 5, originally quotes the scale height as 

H = | [30 + 0.2 (h - 200)] 

3 L m j 

but, based on 292 ionograms at Pennsylvania State University, he reduces this by a 
factor of 1.61 to give equation (6a). 
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It is not known how closely other ionospheres, removed in time or dis- 
tance from Freeman^ sample set, adhere to equation (6a), However, since equation 
(6a) does represent real ionospheres for this region, and according to Wright, a 

similar relationship is valid for other regions, it was used to generate values of 

4 

scale height, H, for these comparisons. The factor of -g in the DC and DODS cor- 
rections for AE, AR, and ARin Tables 2, 3, and 4 is due to our use of equation (6a) 
rather than (6b) in the definition of scale height, H. 

Values of N and h were obtained from each of the midnight and noon 
em m 

mean N g (h) profiles given in reference 11, pages 128 and 129, for Newfoundland, 

Grand Bahama, and Huancayo over a period of 1 year (1959 to 1960) near a solar 
cycle maximum. The lowest, average, and highest of the N em and h^ values ob- 
tained in this manner were then used with equations (5) and (6a) , with a frequency of 
2 GHz, to obtain the low, average, and high values of N^ m and H which are listed 
in Table 5. 

The refraction correction is proportional to N. in all of the formulas 

im 

(except NAP-3 for AE), and hence is inversely proportional to the square of the 
operating frequency. By scaling, the results in this document can be applied at any 
operating frequency; in particular, the scaling factors given in Table 6 should be 
used for the GEOS 1 and 2 tracking frequencies. 

The refraction correction is also proportional to H in the DC, DODS, 
Freeman, and GEOVAP formulas. The scale height, H, which arises from the use of 
the Chapman model in the derivation of these formulas, is determined by means of a 
semiempirical function of the maximum refractivity height,* h . The range of H from 
Table 5 is from 76,667 meters to 150.000 meters; that is, the maximum is about 
twice the minimum value. N. varies by a factor of 10, which means that the product 
R m H can vary by a factor of 20, as shown in Table 5, It follows that for a given ele- 
vation angle, the values of AE, AR, and AR can vary by a factor of about 20 due to 
possible variations in the ionospheric parameters. 

CORRECTIONS FOR AE, AR, AND AR FROM REEK RAY- TRACES 

The REEK ray-trace program, for a given elevation angle, requires as 
input the true range to the target and a refractivity profile. The true range to the 
satellite was obtained for each of the specified true elevation angles from the simpli- 
fied orbit program described in Appendix B, The refractivity profiles used as input to 
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the REEK ray-trace program were generated from all nine combinations of the speci- 
fied values of the N im and H listed in Table 5 using the Chapman model and a compati- 
ble Bent model. 


The Chapman profile is defined as follows: 


N. - N. e 
1 lm 


(1 


Z - e -z ) 



(?) 

( 8 ) 


h. = height of each profile point going from 112.5 km to 1325 km in 
1 12.5 km steps 

The refractivity was set to zero from h^ = 0 to h^ = 112 km. 

The three Chapman and compatible Bent refractivity profiles generated 
from the low, average, and high pairs of N. and H, listed in Table 5 and shown in 
Figure 1, illustrate the shape and possible variation in these profiles. 

The ray-trace values for AE , AR, and AR from the Chapman model are 

listed in Table 7 for all nine of the N. , H combinations given in Table 5. These 

lm 

values are plotted in Figures 2,3, and 4. Several features shown in these graphs are 
of interest and are as follows: 


a. They verify the prior conclusion that the observed variation in N^ m 
has more effect on AE, AR, and AR than the observed variation in H, 
as inferred from equation (6a) and the observed variation in h . 
Because of this , only the low , average , and high combinations of N. m 
and H giveh in Table 5 are used in the subsequent analytic and ray- 
trace refraction correction graphs in Figures 5 through 13. 

b. The ratio of maximum to minimum refraction correction due to a 
change in elevation angle for a given ionosphere is about 100:1 for 
AE , only about 3:1 for AR, and about 10:1 for AR. These generali- 
zations exclude the low probability elevation angles within a degree 
or so of zenith where Ae and Ar both approach zero. Normally AR 
approaches zero at the point of closest approach where E = 0, but for 
an overhead pass (used here and described in appendix B) this occurs 
at zenith. For the chosen 2-GHz frequency and for our average iono- 
sphere, the refraction corrections fall between 0.7 and 0.007 milli- 
degrees (2.5 and 0.025 arc-seconds) for AE, between 9 and 3 meters 
for AR, and between 1.3 and 0.13 cm/sec for Ar. 
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c. At an elevation of about 10°, for both AE and AR, the correction 

curves cross over so that the thicker ionospheres with larger values 

of scale height, H, require a larger refraction correction for &E 

and AR above 10° and less correction below 10° than do thinner 

ionospheres with the same N. but smaller H. 

im 

Table 5 

Low, Average, and High N^ m , h , and H Values (Near Solar Cycle Maximum) 



N 

em 

(10 17 el/m 3 ) 

-N. 

im 

(10“ 6 ) 

h 

m 

(km) 

H 

(km) 

-N. H 
im 

(meters) 

Low (3,3) 

2.19 

2.21 

280 

76.667 

0.169 

Average (2,2) 

10.60 

10.67 

364 

104.667 

1. 117 

High (1,1) 

23.99 

24.14 

500 

150.000 

3.621 


Table 6 

Correction Factors for Other Tracking Systems 


System 

Frequency 

Multiplication 

Factor for N. 

im 


Up (MHz) 

Down (MHz) 

Up/Down 
Equivalent (MHz) 


C-band 

5690 


5765 

5727 

0.122 

TRANET 




108 


342. 936 





(NSSDC data) 






162 


152.416 





324 

> * 

38. 104 





972 ] 


4.234 

Secor 

420.9 


224.5 

280.18 

50.972 


420.9 



434.27 

21.210 




[** 224.5 1 



21.210 


y ■ 




1. GRARR* 







Channel A 

2271.9328 



1928.6 

1.075 

Channel C 

2270. 1328 



1928.0 

1.076 

2. GRARR** 

i 






Channel A 

1799*2 


2253 

1988.3 

1.012 

Channel C 

1801.0 


2253 

1989.5 

1.011 

3. GRARR 







(GEOS-C) 

2074. 6375 


2253 

2158.3 

0.859 

Proposed 







*For GEOS-1 and GEOS-2 
**Between GEOS-2 and GEOS- 

-c 
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Table 7 

REEK Ray-Trace Values From Chapman Model 
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25.247 

21.412 

16.185 
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11.034 

9.805 

AR (em/sec) 

0.4207 

3.1884 

3.7078 

3.1378 

2.2258 

0.0017 

150 

10.67 

? 

w3 
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0.6089 

0.6263 

0.5311 
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0.0000 

g 
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0.1360 

0.1954 

1.4098 

1.6390 

1.3867 
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0.0011 

150 

CJ 

i 

Cd 3 
0=3 

1 

0.1257 

0.1293 

0.1097 

0.0653 

0.0390 

0.0226 

0.0000 

AR (ra) 

2.305 

2.301 

1.951 

1.475 

1.175 

1.006 

0.894 

? 

to 

V. 

•a 

o 

0.0145 

0.0395 

0.2905 

0.3378 

0.2859 

0.2027 

0.0001 

fr- 

ee 

CO 

o 

N 

1 
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I 
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1.5723 
1.1446 
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0.3291 
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0.0000 

? 

« 
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20.411 
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16.419 

11.803 

9.181 
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AR (cm/sec) 

0.0080 
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3.0300 
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2.4186 
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0.0019 

s 

o 
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AE 
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0.0000 

g 
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0.0300 
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0.0000 

g 

03 
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1.076 

0.837 

0.709 

0.626 

AR (cm/ sec) 

0.0154 

0.0420 

0.2760 

0.2786 

0.2203 

0.1514 

0.0001 

76.667 

24.14 
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<D 
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2.21 
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d. The AR corrections for the ionosphere are maximum for an over- 
head pass at elevation angles between 15°, for the thinner iono- 
spheres (low H), and 30° for the thicker ionospheres (high H). The 

G 

AR values decrease for both lower and higher elevation angles in 
sharp contrast to their behavior in the troposphere. The decrease 
in AR for lower elevation angles contradicts the results by Millman 
(reference 32). Weisbrod (reference 33) predicts that the AE iono- 
spheric refraction correction also decreases for lower elevation 
angles. This is also supported by our curves. 

DERIVATION OF COMPATIBLE BENT PROFILES FROM THE CHAPMAN PROFILES 

The Bent profile (reference 23) is given in Appendix A, page A- 10 as 
an exponential function of the parameters and k for the topside ionosphere and as 
a biparabolic function of the parameters N.^ and y^ for the F-layer and bottomside 
ionosphere. The NAP-3 formulas assume the Bent model and therefore incorporate 
the Bent model parameters. The NAP-3 formulas were compared with REEK ray- 
traces through some selected Bent ionospheric profiles as well as the Chapman pro- 
files to determine whether the use of the different models has a significant effect. 

The different model profiles were made compatible by constraining the total electron 
content or total refractivity integral or zenith A R corrections to be the same. 

\ h 

Thus ^ N. (Bent)dh= ^ ^N. (Chapman) dh « N* m He 1 

Since the total electron content is the dominant feature of the ionosphere 
with regard to refraction corrections, this constraint reduces the differences between 
the models to second order effects due to shape differences. 

It is also assumed that the parameters N. and h are the same for both 

* lm m 

models . 

The Bent model topside profile is represented by an exponential decay 
function with scaling coefficient ^ = H^. The Chapman model profile rapidly 
approaches an exponential decay function with scaling coefficient H at heights of 
several H above h^. Therefore, we somewhat arbitrarily chose - H. 

Due to the above assumptions, we have defined the parameters 
Beat [ N im> h m> H b] = Cha P man [ N im’ h m* H ] * 
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The only remaining undefined parameter is the thickness parameter for the Bent 
F-layer and bottomside model. Values of y m were chosen to satisfy the total electron 
content equality constraint above. The effects of this method for defining the compati- 
ble Bent refractivity profiles are shown in Figure 1, where a slight downward shift of 
the compatible Bent profile, relative to the reference Chapman profile, is observed. 

The values for AE , AR, and AR obtained by REEK ray-traces through the 
Bent profiles are plotted in Figures 5 through 13 for comparison with the NAP-3 for- 
mulas and with the ray-trace through the Chapman profiles. 

COMPARISON OF ANALYTIC FORMULAS AND REEK RAY-TRACES 

Values AE , AR, and AR (calculated from the analytic formulas given in 

Tables 2,3, and 4) are plotted in Figures 5 through 13. The values of the Chapman 

ionospheric parameters Nj , h , and H given in Table 5 are used in the DC, DODS, 

Freeman, and GEOVAP formulas. For GPRO the same N. and h parameter values 

lm m 

are used, but the thickness parameter for the GPRO parabolic ionosphere is kept con- 
stant at H =50 km, as in the normal usage of GPRO. For NAP-3, the same N. , 

g ’ ^ ’ lm ’ 

h . and H = H, parameter values are used and, in addition, y is determined as de- 
scribed above so that AR (E = 90°), calculated from Bent’s analytic formula (A27 or 

A28), is the same AR (E = 90°) determined from the REEK ray-traces. The same is 

• o 

true for AR (E = 90 ). The parameter values used are given below: 



-N. 

lm 

(10 -6 ) 

h 

m 

(km) 

H 

(km) 



y m 

(km) 

Low 

2.21 

280 

76.667 

50 

76.667 

172.833 

Average 

10.67 

364 

104.667 

50 

104.667 

235.944 

. High 

24.14 

500 

150.000 

50 

150.000 

337.603 


These same parameters are used in generating the Chapman and Bent pro- 
files for input into the REEK ray-tracing program. The results of these ray- traces 
are also shown in Figures 5 through 13. There is only a slight difference between 
the REEK ray-trace through the Chapman and compatible Bent profiles at low ele- 
vation angles, as can be seen from Figures 5 through 13. 

The DC/DODS correction values are all 20% below the zero order cor- 

• o 

rection AE. , AR. , AR. for elevation angles above 10 , and about 20% below the 
to to io 4 

ray-trace corrections for angles below 10. This is due to the - factor explained 

o 
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earlier in connection with the semiempirical equation for H. If equation (6b) rather 

than (6a) had been used to define H for the Chapman profiles in the REEK ray-traces, 

the DC/DODS corrections for E <10° would have closely agreed with these ray-traces. 

Furthermore, if the DC/DODS formula for E < 10° had also been used for E > 10°, 

the DC/DODS corrections for E > 10° would have better agreed with the ray-traces. 

This can be seen from the dotted extensions of the DC (E < 10°) graphs. Multiplication 

5 

of the DC/DODS curves by the factor g is equivalent to assuming the same total 
electron content for the DC/DODS and REEK ray-trace ionospheres. This was done 
for the selected elevation angles E = 1.5°, 6°, 15°, and 90° to show the resulting 
good agreement between the DC/DODS and the ray-trace results in Figure 5-13. 

The Freeman correction for AR deteriorates rapidly for E < 30°, as 
Freeman acknowledges in reference 5. His correction for A R deteriorates for 
E < 50°. 


The GPRO results are quite good for our low values of the ionospheric 
parameters N. and h , but, for our average and high values of N. and h , the 
GPRO results are consistently too low. This is most likely due to the use of a con- 
stant value of H =50 km for the parabolic model thickness parameter. If the H 
& 5 

parameter were allowed to adjust, to keep the total electron content in the Willmann 

GPRO model equivalent to. that in the Chapman profile used in RE'EK, the GPRO 
results would match the ray-trace results* more closely for all cases . This can be 
accomplished by setting: 


then 


4N. H = e N. H 
im g im 


H = (y) H 
g '4 1 


rather than 50 km 


For the chosen values of H, the new values of H are calculated below 

g $ 

along with the proper multiplicative factors to adjust Ar and Ar. 



H 

(km) 

H = (e/4) H 
g (km) 

Multiplicative Factor for 

AR and AR 
H 

m = 

50 

Low 

76. 667 

52.101 

1.0420 

Average 

104.667 

71. 129 

1.4226 

High 

150.000 

101. 936 

2.0387 
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Applying these multiplicative factors to the GPRO graphs for A R and A R 
for the selected values of E = 1.5°, 6°, 15°, and 90° yields the improved values indi- 
cated in Figures 8 through 13. 

The NAP-3 values for AE 'are consistently lower than the REEK ray-traces 
for all elevations. The NAP-3 values for AR and AR were forced to agree with ray- 
trace values at E = 90°, but are consistently higher than the ray-traces for all other 
elevation angles. Although the magnitudes of the NAP-3 corrections do not agree with 
the REEK ray-traces, the NAP-3 correction curves, in general, match the shapes of 
the ray-trace curves. A modified version of the NAP-3 ionospheric refraction for- 
mulas will soon be available and may improve these results. 

Assuming the same total electron content, it can be shown that N im He 1 = 

40.3 Kp/f 2 , and the only difference between NAP-3 (equation A27) and DC (E < 10°) 

(equation A2a) formulas for AR is in the use of R^ = h - y^ in NAP-3 compared to 

h alone in DC. In Figure A-3, if the line from the center of the earth were drawn to 
m 

intersect the ray-line at h within the ionosphere rather than at h - y at the bottom 
of the ionosphere, the resulting NAP-3 formula for Ar would be identical to the 
DC (E < 10°) formula for AR and would match the REEK ray -trace better, assuming 
the NAP-3, DC/DODS, and REEK ionosphere all have the same total electron content. 
The same conclusions apply to the NAP-3 formulas for AR. * 

We have also tested the preliminary versions of the GTDS formulas . 
However, since these formulas are still under development, the results are not 
reported here. 

DEVIATION AMONG REFRACTION CORRECTIONS 

For elevation angles above about 10°, where most satellite tracking is 
done, the refraction error curves in Figures 5 through 13, deviate in general, from 
each other less than at lower angles . 

The maximum deviations among the corrections for the average iono- 
spheres in Figures 6,9, and 12 are given in Table 8 for selected angles. The Freeman 
corrections are not included below E = 30°, since Freeman specifically warns against 
it in reference 5. As explained earlier, some of the spread among the computed re- 
fraction corrections could be reduced by using different values of the thickness param- 
eters H, H g , H^, y m or by a different choice of effective refraction height of the iono- 
sphere. 
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Table 8 


Maximum Deviation Among Refraction Corrections 


E 

(degrees) 

6AE 

(millidegrees) 

6AR 

(meters) 

6AR 

(cm/sec) 

6 

0.828 

22.34 

29.79 

10 

0.390 

11.30 

11,29 

15 

0.223 

6.27 

5.16 

20 

0.168 

4.09 

2.94 

30 

0.128 

2.30 

1.95 

40 

0.103 

1.60 

0.79 

60 

0.058 

T. 07 

0.35 

90 

0.000 

0.90 

0.00 


RELATIONSHIPS BETWEEN REFRACTION CORRECTIONS FOR AE, AR, AND AR 

NOMINAL ZERO ORDER CORRECTIONS 

Several relations connecting the nominal refraction correction formulas 
for AE, AR, and AR are apparent by inspection of these formulas. They are 

-N. H e 1 AR. 

AE. = — ^ ctn E = - 7 ~ cos E 

io h t h t 

AR. = -N. H e* esc E = h, A E. sec E 
10 im t io 

AR. - +N. He^ E esc E ctn E = - AR. E ctn E = 
io im io 

. • d < AR io> • d < AR io > 

t io dt dE 

HIGHER, ORDER CORRECTIONS 

It is of interest to determine whether relationships similar to the above 
also hold for the higher order GEOS ionospheric refraction correction formuals for Ae , 
AR, and Ar in Tables 2, 3, and 4. The only formula for A E which is more com- 
plicated than the nominal formula is the NAP-3 formula. There seems to be no simple 
relation betweenAE and AR for NAP-3. 
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It can easily be shown that for all the GEOS ionospheric refraction correc- 
tion formulas in Tables 3 and 4, the range rate refraction correction formula for AR, 
is simply the time derivative of the formula for AR, as might be expected for consis- 
tent levels of approximation for AR and AR. 


Thus, 




holds for all formulas. This is not generally stated in the earlier source documents 
where the formulas for AR are derived more laboriously from basic principles. 

The GEOS range rate measurements are all derived by determining the 
number of doppler cycles in a short time interval. The number of doppler cycles in 
this time interval corresponds to the range shift in satellite position relative to 
the station during the same time interval, where range shift is measured in wave- 
lengths of the transmitter frequency. The measured range shift includes the dif- 
ferential effect of the atmosphere on the doppler cycles at the beginning and end 
points of spacecraft position corresponding to the beginning and end points of the 
shift measurement time interval. Therefore, it is reasonable that the range rate re- 
fraction correction should be expressed as the difference in range refraction cor- 
rection divided by the corresponding difference in time, or as a derivative, with 
respect to time, of the range refraction correction formula, 

SUMMARY 

The purpose of this study was to compare the different ionospheric re- ' 
fraction correction techniques used with GEOS E, R, and R data to determine whether 
significant deviations exist among these techniques. If deviations exist, it is difficult 
to know which technique gives the best results, since no perfect model of the iono- 
sphere is available, However, REEK ray-trace results through typical Chapman and 
Bent model ionosphere profiles were included as references for comparison of the 
results from the various formulas and should indicate which formulas are best. 

The following statements summarize the results of this study: 

1. Some of the GEOS ionospheric refraction correction formulas and the 

REEK ray-trace results for AE, AR, and AR are in poor agreement 

for E <30°. The zero order corrections AR. and AR. tend to over- 

10 10 

correct at all elevation angles, becoming worse at the lower elevation 
angles. At E = 10°, AR. q is twice the REEK ray- trace value and 
AR q is ten times the REEK value, so that no correction at all is 
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preferable to the zero order corrections. AE io also tends to over- 
correct for E < 10°, but under-corrects for E > 10° and matches REEK 
at about 10°. 

2. For E >10°, where most satellite tracking is done, the differences in 
AE, AR, and AR are insignificant in comparison to most operational 
requirements. However, for the GEOS Observation Systems Inter- 
comparison Investigation (GOSII), these differences (refer to Table 7) 
are larger in size than observed differences in system biases for 
some of the best electronic geodetic tracking systems. Therefore, 
the systematic differences in refraction correction techniques are 
affecting the GOSH results . 

3. For all the GEOS ionospheric refraction correction formulas ex- 
amined in this report, the range rate refraction correction for AR 
is the time derivative of the formula for AR. In the earlier source 
documents, the AR and AR formulas are derived independently. 

4. Assuming the same total electron content, the NAP-3 and DC/DODS 
corrections for AR and AR are identical except for the use-of an 
effective ionosphere height of R^ = h m - y m in NAP-3 compared to 
h^ alone in DC/DODS. This may account for the differences between 
NAP-3 and REEK. 

5. The REEK ray-traces show that the ionospheric correction for AR 
has a strong maximum between 15° and 30° elevation, contrary to the 
tropospheric correction and contrary to some reports in the literature. 

6. All of the ionospheric correction formulas for AE, AR, and AR 

depend on the integral of the refr activity or electron content. In 

all the ionospheric models discussed here (Chapman, Willmann, 

and Bent) this integral depends on the product of N.^ and some 

thickness parameter. Reasonably accurate values of N. and li 

lm m 

based on bottomside ionograms, are available. The thickness param- 
eter H used in the formulas based on the Chapman model is not well 
known. In fact, the DC/DODS semiempirical linear relationship used 
to derive H from h m differs from the Freeman version by a factor 
of 20%. The uncertainties in specifying the thickness parameters can 
theoretically be reduced by means of topside ionograms or Faraday 
rotation measurements, or multiple frequency techniques such as those 
employed by the Secor and TRANET systems. 
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7. The DC/DODS (E < 10°) corrections programmed for AR and AR 
should also be used for E >10°. The currently programmed correc- 
tions for E > 10° are flat earth approximations which lead to sig- 
nificant errors above 10°. 

8. The DC/DODS elevation angle correction has the proper magnitude 
for a single path angle correction, not the dual path interferometer 
correction intended. The interferometer correction is smaller in 
magnitude, and opposite in sign, and may be obtained from any of 
the range rate (AR) correction equations, as explained in ref- 
erence 1. 

9. All of the formulas compared in this report assume the satellite is 
well above the bulk of the ionosphere. This is a good assumption for 
GEOS-1 and 2 whose perigees are both above 1080 km. However, 
many satellites operate within the ionosphere with electrons above 
as well as below the orbit. In these cases, the ionospheric refrac- 
tion correction formulas based on the total electron content will in 
general be in error to the extent that it may be preferable not to make 
corrections. Some programs, such as NAP-3, provide a separate 
formulation for this case. 
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Figure 1. Chapman and Compatible Bent Refractivity Profiles 
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Figure 2 , Elevation Error Due to Ionospheric Refraction 
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Figure 3 . Range Error Due to Ionospheric Refraction 
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Figure 4. Range Rate Error Due to Ionospheric Refraction 
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Figure 5. Elevation Error Due to Ionospheric Refraction 


24 





ELEVATION REFRACTION ERROR (MILLIOEGREES) 



Figure 6. Elevation Error Due to Ionospheric Refraction 
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ELEVATION REFRACTION ERROR (MILLIOEGREES) 



Figure 7. Elevation Error Due to Ionospheric Refraction 
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RADIO RANGE REFRACTION ERROR (METERS) 



Figure 8. Range Error Due to Ionospheric Refraction 
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Figure 9. Range Error Due to Ionospheric Refraction 
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RADIO RANGE REFRACTION ERROR (METERS) 



Figure 10. Range Error Due to Ionospheric Refraction 


29 



RANGE RATE REFRACTION ERROR (CM/SEC) 



Figure 11. Range Rate Error Due to Ionospheric Refraction 
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RANGE RATE REFRACTION ERROR (CM/SEC) 
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APPENDIX A 

REFRACTION CORRECTION FORMULATIONS 


A.l Different groups using data from the GEOS geodetic tracking systems em- 

ploy different refraction correction equations. This appendix lists some of the differ- 
ent refraction correction equations, and describes how they are being used. For com- 
parison of the equations , the original forms and notations are converted into a common 
form and notation. 

A. 2 GSFC ORBIT DIFFERENTIAL CORRECTION PROGRAM (DC) 

The GSFC Orbit Differential Correction Program (DC) was developed by the 
Advanced Orbit Programming Branch at GSFC. It is used for operational orbit up- 
dating and contains optional refraction correction formulas for interferometer direction 
cosines and for range and range rate. It is documented in reference 2, program 
No. F116. 


A. 2.1 


where 


DC MINITRACK DIRECTION COSINE CORRECTION OR ELEVATION 
ANGLE CORRECTION 


The direction cosine corrections are 


i 



and m 

c 


m 

q 


equation (2), F116 


i, m = measured direction cosines 

0 , m = refraction corrected direction cosines 
*c 9 c 


1 + N 


q =■ 


1 + a - 

s 


N-i (h, t) 


ll t 


equation (4) , F116 


A =N = (n - 1) = tropospheric refractivity at the station; a constant 
s s 

for a given station for a given month. For example , 
a typical N g = 0.000313. 

1 + N = 1 
o 


A-l 



Keeping the ionospheric part 


Am = m - m = 
c 


+mNj (h, t) 


are the changes in observed direction cosines due to refraction in the DC program. 

As shown below, these equations are equivalent to a single refraction cor 
rection equation for elevation angle. 

By definition £= cos E sin A m = cos E cos A 

where E = apparent elevation angle of incident ray 
A = apparent azimuth angle of incident ray 

Refraction causes small changes in elevation angle AE, but in azimuth, 

AA = 0. 

Therefore , taking the derivatives 

A/ = - AE sin E sin A Am = - AE sin E cos A, 


AE 


- iN (h.t) 

-Ai _ 1 ~HT 

sin Esin A sin E sin A 


-m N (h, t) 

AE - - 

sin E cos A sin E cos A 


N (h, t) 

= r ctn E 

h t 

Nj (h, t) = N d H d e V 


N- (h, t) 

— r ctn E 

h t 

equation (19), F116 


H d- kH -t 


[30 + 0.2 (h - 200|] 




V= 1 -e l 



N d-- N lm >0 


equation (8) , F116 

equation (6a) and (6b) 
of this document 

equation (5) , F116 


equation (6) , F116 


A-2 



4 6 

Average GEOS-2 height, h =-=-x 10 meters, and h = 364 km; then 

o m 


H , = 83.73 km, and u = -11.58 
d 

V = 1 - 9 x 10~ 6 = 0. 99999~1 
e V = 2.718255 = ~e 


N 1 (h, t )=N d H d e v »-N lm (|H)e 1 


N (h, t) 

AE = E - E = ctn E = 

C h t 


N. H e l 




ctn E 


AE=-.f- AE. 

5 10 


(Al) 


Thus, both of the equations (2, F116) lead to a single refraction correction 

equation for the elevation angle. Note that the sign is incorrect, that the correction is 

4 

for a single ray and not an interferometer ray pair as intended, and that the factor 

o 

arises from assuming the equation (6b) rather than equation (6a) relation between scale 
height, H, and the F -layer height, h m . 


A. 2. 2 DC RANGE CORRECTION 

From equation (12), F116, for E L 10°, 

AR 1 = R - R c = C A g cse^ + B g N^ (h, t) csc<£ 2 


From equation (14), F116, for E > 10°, 

AR - R - R - C-A esc E + B N- (h, t) esc E 
2 els si’ 

where R = measured range, and R c = refraction corrected range 

For the ionosphere alone with E 4= 10° 

AR, =BN (h, t) csc<£ 

1 s 1 z 

where 

<t> 2 = cos 1 

R = earth radius = 6,378,166 meters 
e ’ ’ 


e cos E\ 

iR + h , 
e my 


equation (17), F116 



B =1 
s 


csc<f> = 


1 - 


^ R e cos 2 
R e + h m, 


- 1/2 


AR, = -N. H . e 
1 im d 


1- 


/R cosE> 

_e 

R + h 
e my 


- 1/2 


Since H, =*4 — H and AR.- » -N. He esc E 
d 5 io im 

For E 4 10° 


AR, = -4-N. H e 1 
1 5i 


/R e cosE \ 2 

\R + h , 
e my 


~ ljf2 AR. 4- sinE 

IQ P 


1 - 


R cos E 

e 

R + h 
e m } 
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(A2a) 


For E >10 

AR 0 =-N. H, e 1 esc E = + | AR. 

2 im d 5 10 


(A2b) 


A. 2. 3 DC RANGE RATE CORRECTION 

From equation (13) F116, for E < 10 C 


. • ® • 

AR, = R - R , = 
1 cl 


r . ->dN, (h, t) 

L B s Wj dh CSC ^2 

2 • 

-CjA g esc <f>^ ctn cos E sin E 


R. 


-B g N 1 (h, t) esc Vtn^lfrTh 


e ms 


(t) 


E sinE) 


From equation (15) F116, for E >10 


* © • • 
ARg = R - R 


c2 


dN x (h, t) 
dh 


[b s £ (t) • r* (t)] 

\ 

- jCjAg + B g N r (h, t)] E ctnE 


esc E 


A-4 



where 


R = measured range rate 

• 

R c = refraction corrected range rate 

The ionospheric part for E 410° is 

- dN, (h, t) 

ARj = B g [? (t) • r_* (t)J gg esc <f> 2 


-B N (h, t) esc 2 ^ ctn <f> L V - h" ) S 

\ e m / 


E sin E 


(n 1 (h, t)^= im ^ - — = "N im e v e u equation (20) , ^116 


for GEOS height; N. e u S 0 

° lm 
For E 4 10° 


Ar = +N H, e 

1 im d 


1 (?^V')® cosEsinE [ 1 '^rrl _ ) 

\ e m/ L \ 6 m/ 


-,-3/2 


cos 2 E 




• 3 tt 

sin E 


R 


-.-3/2 


1 - 


iR + h 

V e mV 


cos E 


(A3a) 


For E >10° 

• • 1 • 

Ar 2 = _n^ (h, t) E ctn E esc E = N. m H d e E ctn E esc E 

= + 4“ AR. (A3b) 
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A. 3 FREEMAN CORRECTION PROGRAMS 

The following formulas were formulated by J. J. Freeman Associates, 
Inc. , under Contract NAS5-9782 for the Operations Evaluation Branch (OEB) at 
Goddard and are documented in reference 5. The purpose of the contract was to de- 
velop better refraction correction equations for use in the calibration of angle, range, 
and range rate systems. 
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A. 3.1 


FREEMAN ELEVATION ANGLE CORRECTION 


The Freeman elevation angle correction (flat-earth approximation) is 


where 


A. 3.2 


' d \) ' hl t 

E = observed elevation angle 

E c = refraction correction elevation angle 

N =(u-l) = tropospheric refractivity at the station 
s 

Keeping the ionospheric part 

ae = - y ht N _ dh 

t O 

Ae = -N^JIe 1 ctn E = Ae. q (A4) 

h t 10 

FREEMAN RANGE CORRECTION 

The Freeman range correction (quasi -flat-earth approximation) is 


/ 


N.dh 

i 


(A 13) 
ref. 5 


AE = E 


E = ctn E 
c 




where 


AR = R 


R c cosff 


2 

J |n| dh - ta ^ ' J j N | h'dh' 


(A25) 
ref. 5 


ip o 

R = electromagnetically determined range 

R c = true refraction corrected range (meters) 

h = satellite height ss 10^ km for GEOS 

h’ = distance from earth surface to ray path 

R = distance from earth center to surface 
e 

cc= ( n / 2 - E) = apparent zenith angle 
I N | = absolute value of tropospheric or ionospheric refractivity = |u-l | 

h oh 


AR = 


sinE 


/ M 


dh 


ctn^E 

R_ 


/ i N i 


hdh 


In the previous equation, if |n| assumes the value -N^ 
then the ionospheric formula is 


Ar = -N. H e 

lm 




esc E 


(A5) 
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A. 3.3 


FREEMAN RANGE RATE CORRECTION 


From equation (A51), reference 5 



where 


V^RE 


As h T approaches the GEOS satellite height of approximately 10 km, the 
refractivity at the target, N^, approaches zero. Therefore the range rate correction 
equation reduces to 


AR =- 


E cosE 
sin 2 E 


/ 1 


N I dh 


,( sinE ] j 1 |n| h dh 


R 


In the previous equation, if I N assumes the value -N 

i 

f H + h 


AR =+N. H e 
1m 


4 -— 3 —) ('■ 

sin E [ \ sin E / \ 


m 


R^ 


-) 
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A. 5 GSFC GPRO PROGRAM 

The GEOS Preprocessor Program (GPRO) was developed for the OEB at 
GSFC by DBA Systems, Inc. , under contract NAS 5-9860, and is documented in ref- 
erences 4 and 22. It is used with the Geodetic Data Adjustment Program (GDAP) , and 
with the Network Adjustment Programs (NAP-1 and NAP-2). The GPRO program 
utilizes a parabolic model for the ionosphere. 



A. 5.1 


GPEO RANGE REFRACTION CORRECTION 


The GPRO range refraction correction, from reference 4, appendix C, 
section la, is 

20 ! 2 h 

T 5 T . r _._2 „ , in 2 2 P 11/2 


AR, = 


^ 1 sin E^ + sin^Ej + 4/?^ cos^Ej 


2j3 1 esc E 1 


1 + ^1 + 4/3 2 2 ctn 2 eJ 1/2 

where AR. = correction for ionospheric refraction 


E^ = local elevation angle of ray upon entry to ionosphere 


The angle E^ can be computed from 


cos E„ = 


R cos E 
e 


1 R + h - 3 H 
e m g 


where 


E = apparent elevation angle at observer 
h m -3Hg= height of the base of the ionosphere 


/Sg are given by the following functions of ionospheric parameters 


/}, = -2H 

P 1 . g 




80.5 N 


m 


-= - 4 H N. 

g im 


o 25 H 

4 fT — . ■ .... — £ 

p 2 3 (R + h - 3H ) 
e m g 


where 

f = ranging frequency (or modulation frequencies) 

f = critical frequency (f 2 = 80.5 N ) 
m m m' 

Nni = maximum electron density 

XT f 1 /f m \ 2 -40.25 N im 

N im = maximum refractivity, N im = p — — 

h m = height corresponding to maximum electron density 

Hg = thickness factor of parabolic ionosphere = 50 km 

R = radius of the earth 
e 
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ARj =' 


-8 H N. esc E. 

K , un 1 


1 + 


25 H 

1 + F(r J 

L ' e 


ctn 2 E, 

J 5 i 


TT /2 


+ h - 3 H 
m 
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A. 5. 2 GPRO RANGE RATE CORRECTION 

The GPRO range rate refraction correction, from reference 22 is. 


where 


A . - 2 ^D 

AR = 

D, 


2 2 2 
D^= sin + (sin E^ + 4/3g cos E.^) 


1/2 


- 1/2 


• r 2 2 2 ' 

D.j= ccfe E^ + (sin E^+4^ cos E^) (sin E.^ cos E^- 
4^2 cos E 1 sin E^J E^ 


R cos E 
Cos E i = R + h 

m • I 


l e m 
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(R g sin E) E 


1 (R + h - 3 H ) sin E 
'em g' 1 


- 2 / 3 , 
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cos E^ + 


sin E^ cos E^ 


4/^ cos Ej^ sin E^ 


. . 2 ' . . n 2 2 „ , l/2 

(sin E^ + 4/?2 cos E^) 


2 221/2 
sin E.^ + (sin E.^ + 4 ( 3 ^ cos E^) ' 


AR 


• 
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(1 + 4/J 2 ctn 2 E.) 1 ^ 2 
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-8 N 


lm 


H E, ctn E, esc E, 
K 1 1 1 


AR =- 


1 + 


(1 + 4^ ctn 2 E x ) 1/2 


1 + ( 1 + 4/? 2 ctn 2 E x ) 1/2 


A. 6 GSFC NAP-3 PROGRAM 

The Network Adjustment Program (NAP-3) formulas were formulated by 
DBA System, Inc, (Contract No. NAS5-17730), and are documented in reference 23. 
Further explanation of the Bent profile, NAP-3 Elevation Angle Correction, NAP-3 
Range Correction, and NAP-3 Range Rate Correction is presented in the following 
paragraphs (reference 23 and supplementary notes by Dr. R. Bent). 

A. 6. 1 BENT PROFILE 

The Bent profile was obtained from the analysis of many tens-of-thousands 
of topside and bottomside profiles and it fits closely to actual profiles. It contains a 
biparabolic lower portion and an exponential upper portion and is used in all NAP-3 
computations on satellite refraction effects. The Bent profile is illustrated in Fig- 
ure A-l. 



where 


Figure A-l. Bent Profile 


k = exponential decay constant 
N = electron density at height h 


N = maximum electron density 

y = thickness factor for biparabolic 
m model ionosphere 


A. 6. 2 


NAP-3 ELEVATION ANGLE CORRECTION 


The passage of a ray path through the ionosphere is illustrated in Fig- 
ure A-2 from whence it can he shown that the deviation angle a = = , say, a . 

The main angles under consideration are the angle of elevation E^ the de- 
viation angle a ^ (or a ) and the error in the elevation angle p . The angle at the sat- 
ellite between the line of sight and actual path is (j3- a). 

First, obtain the deviation angle a and from this compute /3. 


The deviation angle - Q^ f but as let them all be 


called a :• 


where 


a = $ T - 8 l 


(A9) 


= the angle subtended at the center of the earth by the limits of 
the true ray path in the ionosphere 


where 


0^ = the angle subtended from the same height limits by the ap- 
parent ray path. 

Snell's law gives: 

p p sin i = p sin <p = constant (A10) 


i and <P = the angles of incidence of the true and apparent ray paths 
respectively at height p , the following equation is derived: 


d 0=ta^l dp 


(All) 


which, when combined with equation (A10) yields the following equation: 


d0 = 


sin <f> 


p J - sin^ <P 


dp 


(A12) 


Simplifying the Appleton-Hartree equation and neglecting magnetic field effects 
(which can safely be done at very high frequencies) the equation for the index of re- 
fraction is given by 


P 


2 


= 1 - 



(A13) 


where f = the frequency of transmission f^ = the plasma frequency at the height 

under investigation. 


A-li 





However, as the electron density N is proportional to f 2 , equation (A 13) can be re- 
written including the critical frequency as: 


where 



(A14) 


f = the critical frequency of the F layer corresponding to the 
c 

electron density 

N = the electron density at the height under investigation 


Rewrite equation (A14) as follows: 

M 2 = 1 - X 


(A15) 


Combining equations (A12) and (A15) yields: 


$ 

where 


-j 

^ a 


sin <P 


4 


dp 


Pj cos - X 


(A 16) 


a and b = the lower and upper heights of the ionosphere. 


To obtain Q ^ of equation (A9), assume that the ray travels in a straight line with 
M = 1 and x = 0. Equation (A16) then yields: 


Q 


l 



tan 

P 


dp - 


(A17) 


Now the frequencies of vertical and oblique incidence radio waves are related by the 
expression: 

f = f sec <p 
v 

where 

f = the frequency of the oblique incidence wave, 

f = the equivalent vertical incidence frequency, and 

= the angle of incidence of the apparent oblique incidence 
ray path. 
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Now at the point of reflection 7 — sec <P = 1 and if this function is less than unity the 

1 fw 

ray proceeds through the ionosphere and therefore the value of sec is a measure 
of the deviation of the penetrating ray. Its maximum value will occur at or just below 
the height of maximum electron density. An increase in this deviation factor indicates 
an increase in the deviation angle of the penetrating and apparent. ray path. The maxi- 
mum deviation occurring when the value of -y- sec is just below unity. Let this de- 
viation factor be u where: 


/ f N f 

u = r j — sec v? J 


which, from equations (A13) and (A15) gives 

u - x sec 2 <P. 

Therefore equation (A16) can be rewritten as: 

• b 


e n 


-j 


tan <P 
1 

P(l-u ) 2 


dp 


hence 


-« T - / b ^[<i-ur7-i]d P . 


(Ai 8 ) 


(A19) 


At the height of refraction u = 1 and at the height of maximum electron density u is 
simply the square of the deviation factor and for all our purposes u < 1 . 

For all rays that penetrate the ionosphere we can apply the binomial theorem to equa- 
tion (A 19), to give: 


2 + If u 3 + * • -l d P 


hence 


1 fk u tan <p f- , 3 ,5 ! 

2 J a P [ 4 u 8 u 

= 1 r b X ten f sec 2 4i + 3 u + | u 2 + || 
2j a p L 4 8 64 


1 


u 3 + 


Applying the second law of the mean for integrals: 


1 t f k x tan <p sec 2 , 

= 2 ^J p d(> 

a 


dp . 


(A 20 ) 
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where 


£ = a number between the greatest and least values of the 
series over the interval (a, b). 

u = 0 at the limits (a, b) and has a maximum value, less than unity, near the height of 
maximum electron density. 

The value of £ must therefore lie between unity and the value of the series at maxi- 
mum electron density and is defined as the mean of these two values. 

Applying the second law of the mean for integrals once more to equation (A20) yields: 

' 2 0 f b 

x dp (A21) 


, tan sec' 

a =- £ 2 

2 S P n 


> 


where 


From Figure A2 


then 


p = P (1 + 0.5333-75 — ) = some radius between limits a 
0 m Pm and b 

p m=% + h m 

<P Q = the value of <P at the radius p Q 
He 


sin = 
o o 


max = 


■ cos E 

(A 


1 - 



e 1 COS 2 E 


(A22) 


(A23) 


£mean is obtained from a table of u max vs ^mean values ‘ ^mean “ 1 £or f = 2 GHz ‘ 


Now, 



(A24) 


Therefore, the deviation angle is given by 


a = 


1 

2 


9 2 

f tan <p sec <P N 

f ) * Po * 


(A25) 


f b 

where N is given by / N d p and is the total electron content of a vertical column 
of unit cross-section of the ionosphere. 
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It can now be shown from Figure A-2 that 


where 


/3 = AE = cos 

X x = (rJ - R 2 cos 2 E) 1/2 + R g cos E tan ) 

X„ = R„ sin E - R cos E tan -5- 
2 e e 2 

R = R q + h 
s e 


X 1 cos CC - X 2 

(X ^ + X g - 2X X X 2 coscc) 1 ^ 2 


(A26) 


R = radius of the earth 
e 

h = height of the satellite 

V 

A. 6. 3 NAP-3 RANGE CORRECTION (REFERENCE 23) 

Assuming the total electron content (N^,) has been computed in the vertical 
direction through the ionosphere along AB (Figure A-3) , then the one-way value of 



5 29- A-3 

Figure A-3. Range Correction Diagram 
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If E is the satellite elevation angle: 


cos 


-HW 


cos^ E, 


where R e is the radius of the earth and is the height of the bottom of the ionosphere 
above the surface, (R^ = h^ - y m ), 

40.3 


AR = 


f 


/ ' ( R e + \ ) 


cos 2 E 


(A27) 


When the satellite height, h , is above the top of the biparabolic layer (h the 

~ h s m & 

equation = / Ndh gives: 

* J o 


n t •/ 
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77T N e da + 
16 m 
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•• N T N ecl 1 16K ( 1 6 I 480 y m 

From reference 23, equations 1, 3, and 4: 
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we have 



(A28) 


A. 6.4 NAP-3 RANGE -RATE CORRECTION 


Equation (A27), when differentiated with respect to time, provides the 
range-rate correction. 


AR = 


Mm 
a e 


e + - 3W 
S h s 





- Mm e 

3 E 


since 


h g = 0 for our circular orbit, 



(A29) 



APPENDIX B 

ELEVATION RATE DETERMINATION 



APPENDIX B 

ELEVATION RATE DETERMINATION 


The required expression for elevation rate is derived from the geometry 
shown in Figure B-l. In this figure a circular overhead orbit around a spherical, sta 
tionary earth is assumed. 


TRACKING STATION 



Figure B-l. Plane Containing Satellite, Tracking Station, and Center of Earth 


t-* R m sin cot 

R cos E t. T 

sin cot = — ; cos E = 


(Bl) 


(R + R sin E) 


cos cot = 


*; sin E = 


(R t cos cut - R g ) 


R 


9 2 l/2 

R = (R + R™ - 2B rr R cos cot) 1 
si is 


(B2) 

(B3) 


substitute (B3) into (Bl). 
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50 



22 - 1 /2 
cos E = R t sin wt (R g + R T - 2R T R g cos wt) 


then by differentiation, 


dE -1 -3 2 

sin E = w R,p R cos wt - w R,j, R (R^, R g sin wt) 


(B4) 


substitute (B2) into (B4) 


-1 dE -1 2 2 “2 

-R (R T cos wt - R g )-gp = wR T R (R cos wt - R T R g sin wt) (R ) 


dE 

dt 


= [- wR t (R 2 coscut + R 2 coscut - 2 R t R g cos 2 cut 


" R T R s + R T R s cos 2 cut) J [r 2 (R t coscut - R g ) 


]" 


dE [- <uR t (R t - R g cos cut) (R T cos cut - R g )] 
dt o — — 


R (R t cos cut - R g ) 


dE 


,2 . 2 


^ “ w Rrp (Rrp * Rg cos wt) (Brp ^ Rg 2Rrp R g cos wt) 
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R t = distance from the center of earth to the satellite (7,711,499 meters) 


R = distance from the center of earth to the station (6,378,166 meters) 
s 


R 


= distance from the station to the satellite 

2 n 


= angular velocity of the satellite = 

R 


T 

T 


= period of the satellite = T 


m 


T 

3/2 


(page 96 of reference 18) 


= 84.347 minutes, the period of a circular’ orbit at the earth's surface 
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